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ABSTRACT 

We present results from 13776 radial-velocity measurements of 1278 candidate members of the old 
(4 Gyr) open cluster M67 (NGC 2682). The measurements are the results of a long-term survey 
that includes data from seven telescopes with observations for some stars spanning over 40 years. For 
narrow-lined stars, radial velocities are measured with precisions ranging from about 0.1 to 0.8 km s“^. 
The combined stellar sample reaches from the brightest giants in the cluster down to about 4 mag¬ 
nitudes below the main-sequence turnoff (U = 16.5), covering a mass range of about 1.34 Mq to 
0.76 M0. Spatially, the sample extends to a radus of 30 arcminutes (7.4 pc in projection at a distant 
of 850 pc or 6-7 core radii). We find M67 to have a mean radial velocity of -1-33.64 km s“^ (with 
an internal precision of ±0.03 km s“^) well separated from the mean velocity of the field. For stars 
with >3 measurements, we derive radial-velocity membership probabilities and identify radial-velocity 
variables, finding 562 cluster members, 142 of which show significant radial-velocity variability. We 
use these cluster members to construct a color-magnitude diagram and identify a rich sample of stars 
that lie far from the standard single star isochrone, including the well-known blue stragglers, sub¬ 
subgiants and yellow giants. These exotic stars have a binary frequency of (at least) 80%, more 
than three times that detected for stars in the remainder of the sample. We confirm that the cluster 
is mass segregated, finding the binaries to be more centrally concentrated than the single stars in 
our sample at the 99.8% confidence level (and at the 98.7% confidence level when only considering 
main-sequence stars). The blue stragglers are centrally concentrated as compared to the solar-type 
main-sequence single stars in the cluster at the 99.7% confidence level. Accounting for measurement 
precision, we derive a radial-velocity dispersion in M67 of 0.80 ± 0.04 km s“^ for our sample of single 
main-sequence stars, subgiants and giants with V < 15.5. When corrected for undetected binaries, 
this sample yields a true radial-velocity dispersion of O.bO^Qgg km s“^. The radial distribution of the 
velocity dispersion is consistent with an isothermal distribution within our stellar sample. Using the 
cluster radial-velocity dispersion, we estimate a virial mass for the cluster of 2100jl|5Q 

open clusters and associations: individual (NGC 2682) - binaries: spectroscopic - methods: obser¬ 
vational - techniques: spectroscopic - blue stragglers - stars: kinematics and dynamics 


1. INTRODUCTION 

M67 (NGC 2682) is one of the few old and rich 
open clusters in our Galaxy, and therefore is central 
to our understanding of stellar evolution, stellar dy¬ 
namics and star cluster evolution. M67 is located at 
a = 8'‘51™23!3, J = ±11°49'02" 1420001. With an age 
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ation time of about 100 Myr (jMathieu & LathamI 
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[TMl), M67 is a highly dynamically evolved system. 
Importantly, M67 is relatively nearby, with recent 
distance measuremen ts ranging f r om between about 
800 pc to 900 pc llaii^ 119^51: iNissen et al.l [r987t 


Montgomery et afTllOOSt ICarryo et al.l 1994 : Fan et alJ 


1996t Grocholski fc Saraiedinil l2003t ISandauistI 

^ l2007t iPasQuini et al.f 


2004 : 


2008 


Saraiedini et al.l l2009h . and has low extinction, with 
recent E{B — V) determi na tions between 0.015 
and 0.056 d .Ianes fc SmithI 11984 iBurstein et alJ 
19861: iMontgo^rv et al.l 119931: ICarraro et al.l 11994 
Fan et al.l Il996t iTavlod 1200711 . Furthermore, most 
studies agree that M67 has a roughly solar metal- 
licity, with recent [Fe/H] value s ranging from -0.10 
to ±0.05 (I Janes fc Smith 11984 fBurstein et al.l Il986t 
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Tavlod 120071: iFriel et al.l 120101 : iPancino et alJ 120101 : 
Jacobson et al.l l201lh . Thus M67 provides a large 


sample of solar-type dwarf and evolved stars that are 
easily accessible to a variety of ground-based (and 
space-based) observations. 

Indeed, M67 has been extensive ly studied through 


photometry from X-ray s to ne a r-IR feello ni et a.1.] 


VandenBerg &: Stetso: 


ays to ne c 

d [200l 


[Landsman et af 


1998 


1998 


Nissen et aLl 1198'^ Montgoinerv et al. 1993 : Fan et alJ 

19961 Balaguer-Nunez et al.l 120071 : lYadav et al.l 120081 : 
Saraiedini et al.l I2009H . includi ng several tim e- series 


optic al photometric sur v eys (I Gilliland et al.l 
Ivan den Berg et ^ 120021 : IStassun et al.l 


Sandauist fc Shetrond l2003bl: [Bruntt et al.l _ 

Pribulla et al.l 120081 : lYakut et al.l 120091) . These surveys 


have revealed many intriguing stars, some of which lie 
far from the standard single-star evolutionary sequence 
in a color-magnitude diagram (CMD). 

In order to interpret these observations, kinematic 
membership probabilities are paramount. Th ere have 
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; iGirard et al.l 119891: iZhao et al.l(l993l: lYadav 
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. and see alsolLoktinll2005h. to various limitinff 
ies and spatial extents. There have also been 
d-velocitv (RV) survevs (IMathieu et al.l 119861. 

mag- 
a few 

1990; 

Milondll992l:IMilone &: Lathamlll994 lYadav et al.1 
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firmed the kinematic cluster memberships of a rich pop¬ 
ulation of blue straggler stars (BSS; residing blueward of 
and generally brighter than the main-sequence turnoff), 
yellow giants (residing between the BSS region and the 
normal giant sequence), and “sub-subgiants” (residing 
to the red of the main sequence but fainter than the 
subgiant branch, and also known as “red stragglers”), 
among the kinematic cluster members (see Figure [5] and 
Section 16.2|) . 

Numerous theoretical efforts have aimed to explain the 
origins of these exotic stars th rough studies of the dy¬ 
namical evolu t ion of M67 (e.g. iLeonard fc Linnelll 119921 : 
iLeonardI 119961 : iHurlev et al.l 1200511 . These studies show 
that close stellar encounters, and particularly those in¬ 
volving binary stars, may be relatively frequent in M67 
and can lead to the creation of exotic stars similar to 
those observed in the true cluster. Furthermore, these 
models emphasize the importance of binaries to the dy¬ 
namical evolution of the cluster. 

To date, the published results of the binary population 
in M67 has been l imited. The la r gest s tudy of binaries 
in M67 is that of iMathieu et al.l (1199011 . who observed 
a sample of bright (V < 12.8) proper-motion members, 
and present orbits for 22 spectroscopic binaries. We have 
continued to monitor these and other stars in M67 in or¬ 
der to extend our sample of detected binaries (and those 
with orbital solutions) to longer orbital periods, fainter 
magnitudes and a larger distance from the cluster center. 
A progress report summarizing the characteristics of 85 
spectroscopic binaries was pre sented at the General As¬ 
sembly of the lAU in Prague (|Lathaml 1200711 . 

Here we present results from our ongoing RV survey of 
the cluster. To date we have obtained 13776 RV mea¬ 
surements^ of 1278 stars in M67 with 8 < V < 16.5 


^ Throughout this paper, when quoting numbers of RV measure¬ 
ments, we provide the number of RVs from the primary stars (and 


(about 1.34 Mq to 0.76 M0) and extending spatially 
to 30 arcmin in radius from the cluster center (7.4 pc 
in projection at a distance of 850 pc, or approximately 
6 to 7 core radii).^ Our stellar sample spans from the 
brightest stars in the cluster down to 4 magnit udes below 
the m ain-sequence turnoff. Also for reference, iFan et ahl 
(|1996ll find a half-mass radius for stars in our observed 
magnitude range in M67 of about 10 to 11 arcmin. Tidal 
radius estimat es for the cluster range from 50 arcmin 


to 100 arcmin (iKeenan et al.l 

19731: iPiskunov et alJl2008l: 

iDavenoort & SandQuistll201C 

). 


In Section [21 we define our stellar sample in detail, and 
in Sections [3] and m we describe our observations and the 
completeness of our data. In short, our observations are 
nearly complete within our “primary sample” of stars 
with V < 15.5 and within 30 arcmin from the cluster 
center. Within this primary sample, we have at least 
one RV measurement for all but one star and at least 3 
RV measurements for all but four stars (two of which are 
rapid rotators). None of these four stars with < 3 RV 
measurements in our primary sample are proper-motion 
members; moreover, we have at least three RV measure¬ 
ments for all proper-motion members in our primary stel¬ 
lar sample. Our time baseline of observations for some 
stars reaches to over 40 years. 

Throughout this paper, we use a cutoff of at least 
3 RV measurements before attempting to derive RV 
membership probabilities and vari ability statistics . Ex- 
tensive Monte Ga r lo an alyses by iMathieuI (Il983l) and 
I Geller fc MathieiJ (j2012h indicate that 3 RV measure¬ 
ments over the course of at least one year is sufficient to 
detect nearly all binaries with orbital periods of less than 
10^ days. With our observations, we can detect binaries 
with significantly longer periods, out to ~ 10^ days. The 
detection of binaries is particularly important for deter¬ 
mining reliable RV membership probabilities. In general, 
binaries without RV orbital solutions do not yield pre¬ 
cise RV membership probabilities, due to their unknown 
center-of-mass ( 7 -) RVs. We discuss our identihcation 
of binaries and cluster members in Section |5l and we 
present a summary table of our results for each observed 
star in Table [H 

Finally, in Section jS] we use our cluster members to 
investigate the GMD (cleaned from field star contami¬ 
nation), a few stars of note, the projected radial distri¬ 
bution of different cluster populations, and the velocity 
dispersion and virial mass of the cluster. Then in Sec¬ 
tion |7| we summarize our results. 

2. STELLAR SAMPLE 

Our RV su rvey of M67 beg an as part of the disserta¬ 
tion work of IMathieuI (|1983ll . ta king advantage of th e 
GfA Digital Speedometers (DS; iLathaml 1198^ Il992fl . 
Three nearly identical instruments were used, initially 
on the MMT and 1.5-m Tillinghast Reflector at the Fred 

therefore count the number of spectra that result in at least one 
RV measurement). We do not add to our count additional RVs 
from, for instance, the secondaries of SB2s. 

^ IZh ao et al .] 11199611 derive a core radius for M67 of 5.2 arcmin. 
IBonatt^^^BiSl ||200^ use 2MASS data to der i ve a c ore radius for 
M67 of 4.86 arcmin. IDavenport &: Sandguistl II2010I) revised this 
result by reanalyzing 2MASS data, and find a core radius of 4.12 
arcmin (but a much larger core radius of 8.24 di 0.60 arcmin when 
using the fainter stars in their SDSS sample). Therefore our sample 
extends to between approximately 6 and 7 core radii. 
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TABLE 1 

Time Span of Data from Each Telescope 


Telescope 

HJDstart 

HJDend 

Days 

N stars 

.^obs 

Palomar Hale 5m 

2440952 

2445297 

4345 

112 

367 

Tillinghast 1.5m + DS 

2444184 

2454958 

10774 

357 

4670 

CORAVEL 

2444340 

2446413 

2073 

5 

92 

MMT -1- DS 

2445337 

2450830 

5493 

376 

1889 

Wyeth 1.5m -|- DS 

2445722 

2453433 

7711 

34 

332 

WIYN -1- Hydra 

2453386 

2456709 

3323 

1163 

5753 

Tillinghast 1.5m + TRES 

2455143 

2456801 

1658 

111 

672 


Lawrence Whipple Observatory on Mount Hopins, Ari¬ 
zona, and then later on the 1.5-m Wyeth Reflector at the 
Oak Ridge Observatory in the Town of Harvard, Mas¬ 
sachusetts. Subsequently the M67 target samples were 
expanded several times. RV measurements from other 
programs were integrated into the database, and our ob¬ 
servational facilities were extended to include Hydra at 
the WIYN Observatory^ and TRES on the Tillinghast 
Reflector. In Table [T] we list the different telescopes and 
Instruments used for this project, along with the dates, 
number of stars observed and number of observations 
from each telescope. For the record, here we briefly re¬ 
view the history of what now constitutes more than 40 
years of observations, and end by defining the M67 stellar 
sample that comprises the foundation of this paper. 

2 .1. CfA 

After exploratory observations of a few stars, t he initial 
CfA s ample was defined in 1982 as all stars with lSandersI 
(|1977f) proper-motion membership probabilities greater 
than 50%, V < 12.8 and {B — V) > 0.40. This sam¬ 
ple comprised the top of th e main sequen ce, subgiants 
and the giant branch. (The LSandersIIl977l study covers 
a square of approximately 80 arcmin on a side, centered 
on the cluster.) A study to monitor the RVs of the 13 
classical blue stragglers in M67 was being pursued in¬ 
dependently (iLatham fc Miloii3ll996f) . and the relevant 
results from that survey are included In this paper as 
a convenience to the reader. Also, a subset of 28 gi¬ 
ants and subgiants were observed more intensively from 
Spring 1982 to Spring 1983. 

By Spring 1988 t he targe t samp le expanded substan¬ 
tially to include all ISandera (|1977l l members to V = 14, 
and all members in the cluster core (radius < 10 arcmin) 
to y = 16 (only observable with the MMT), totaling 
432 stars. The last surviving CfA Digital Speedometer, 
on the 1.5-m Tillinghast Reflector, was retired in the 
summer of 2009. Over the following five observing sea¬ 
sons the new Tillinghast Reflector Echelle Spectrograph 
(TRES) was used to continue the RV observations of tar¬ 
gets (mostly binaries) from both the CfA and the WIYN 
samples. Whenever the observations accumulating for a 
target suggested that the velocity was not constant, ad¬ 
ditional observations were scheduled at a frequency de¬ 
signed to reveal its orbital parameters. For some systems 
the CfA observations span almost 35 years. 

Importantly, Roger Griffin and James Gunn also had 
a RV program for M67 from 1971 to 1982 at the 

® The WIYN Observatory is a joint facility of the University 
of Wisconsin-Madison, Indiana University, the National Optical 
Astronomy Observatory and the University of Missouri. 


Palomar Hale 5m telescope, which they supplemented 
with contemporaneous observations obtained with the 
CORAVEL Instrument at Haute Provence for five of 
the binaries. Their target sample was very similar to 
our initial sample, and the combination of the two data 
sets to expand the time baseline was a natural and 
straightforward step. The Integration of the Palomar 
and CO RAVEL data w i th th e CfA data is discussed in 
detail in iMathieu et al.l (1198611 , where the entirety of the 
Palomar data as well as the CfA data to that date are 
presented, and in iMathleu et al.l (|199(1tl . 

Hereinafter, we will refer to the target stars and the 
measurements taken at all telescope/lnstrument pairs ex¬ 
cept WIYN/Hydra as the CfA sample and data. As of 
February 2015, there are 447 stars In the CfA sample. 
The radius - V-magnltude distribution of the stars ob¬ 
served by these telescopes is shown in Figure[T](top), and 
the middle and bottom panels of the same figure present 
the chronology of observations graphically. 

2 .2. WIYN Stellar Sample 

As evident in Figure [U the CfA sample Is not com¬ 
prehensive in radi us for th e faint er stars in the sample. 
Furthermore, the ISandersI (jl977tl proper-motion study 
was progressively more Incomplete and less precise with 
increasing magnitude. The Hydra Multi-Object Spectro¬ 
graph (MOS) on the WIYN 3.5m telescope, combined 
with modern target lists, has been able to provide a 
complete magnitude and spatially limited sample, as de¬ 
scribed here. 

WIYN observations of M67 began on January 15, 
2005 as part o f the WIYN Open Cluster Study (WOCS; 
lMathieull2000ll . The current WIYN target list contains 
all stars with 10 < V < 16.5 and within 30 arcminutes in 
radius from the cluster center. We drew our sample (In¬ 
cluding astrometric positions) from t he 2MASS catalog, 
supple mented with photometry from [Montgomery et al.l 
(|1993ll . The Hydra MOS on WIYN has a 1° field of view, 
which sets the radial extent of the WIYN survey. 

We can derive reliable RVs with our WIYN observing 
setup for stars with {B — V)o > 0.4. (Rapid rotation and 
a diminished number of absorption features often hinders 
our ability to derive precise RVs for earlier-type stars.) 
The only stars In M67 that are bluer than this limit 
and bright enough to be within our sample are BSS. Be¬ 
cause of their scientific interest, we included these stars 
in the WIYN stellar sample. In total there are 1278 stars 
within the WIYN stellar sample; 382 of these stars are 
also In the CfA stellar sample. We show the radius - 
V-magnltude distribution for stars observed at WIYN In 
the top panel of Figured 
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Fig. 1.— The CfA stellar sample, showing (top) V magni¬ 
tude as a function of radius from the cluster center, and the 
chronology of the sample, with V magnitude (middle) and 
radial distance from the cluster center (bottom) plotted as 
functions of the HJD of the first observation, respectively, for 
all CfA observed stars. 

As seen in Figure[2] (middle), the V magnitude range of 
the stars observed at WIYN has evolved with time. Ini¬ 
tially we chose to make use of WIYN’s advantage at faint 
magnitudes and prioritized stars from 12.5 < V < 16.5. 
(The Hydra MOS has a dynamic range limit of 4 mag.) 
Later, to improve the completeness of our observations 
within the combined WIYN and CfA sample, we ex¬ 
tended our WIYN sample to include all stars within a 
30 arcminute radius from the cluster center and with 






r 


• 

1 4 • 


r 


• • 

:»' ■ 

t 

i‘- • 

* j 
j ^ 

1 

1 

- 

»: 

< • 

i • 

- 


•* 

•1 

f ■ 

il * 



•• 

J - 

1 



•1 

•• 

* *• 

f * • ! 



t, 
• • 

<: • • 

• 

• 

• 

• 

• 


.il 

1, . 




5.3 5.4 5.5 5.6 5.7 

(HJD-2400000)o/10" 


Fig. 2.— The WIYN stellar sample, showing (top) V magni¬ 
tude as a function of radius from the cluster center, and the 
chronology of the sample, with V magnitude (middle) and 
radial distance from the cluster center (bottom) plotted as 
functions of the HJD of the first observation, respectively, for 
all WIYN observed stars. 


10 < H < 16.5. For a few epochs, we focused exclusively 
on the brighter stars in our list in order to build this 
complete sample. 

We also note that we removed from our sample the few 
high-proper-motion stars in the M67 field whose coordi¬ 
nates have changed significantly during the course of our 
survey, as these stars are certainly not cluster members. 
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3. OBSERVATIONS 

Details about the telescopes, observing procedures, 
and data reductions of spectra obtaine d with the CfA 
Digit al Speedometers can be found in iLathaml (I1985L 
Il992fl . The corresponding information for spectra ob¬ 
tained with Hydra at the WIY N O bservatory can b e 
found in iGeller et al. I ((20081120101) and iHole et al.l (|2009D . 
We describe briefly here the TRES instrument and ob¬ 
servations. 

TRES is a stabilized fiber-fed echelle spectrograph with 
a CCD detector and resolution of 44,000. Procedures 
were adopted to ensure that the RVs from TRES could 
be adjusted to the native velocity system of the CfA Dig¬ 
ital Speedometers. Although TRES delivers wavelength 
coverage from 390 to 900 nm, only the order centered 
on the Mg b features was used to derive RVs in order 
to match the spectral region used by the CfA Digital 
Speedometers, using the same library of synthetic tem¬ 
plate spectra for the correlation analysis. To establish 
the zero point offset between the TRES and CfA Digi¬ 
tal Speedometer RVs, observations of the afternoon blue 
sky and of selected lAU standards were obtained on most 
nights. There were several modifications to TRES during 
the first years of its operation, and corresponding zero- 
point shifts as large as 0.1 km s“^ were measured. Since 
March 2012 the TRES zeropoint has been stable at the 
level of about 0.01 km s“^. 

The TRES velocities are all (initially) shifted by the 
gravitational redshift of the Sun and blueshift of the 
Earth, because the library of templates does not in¬ 
clude either of those effects. Therefore the derived stellar 
velocities are all redshifted by the net amount, nomi¬ 
nally -1-0.62 km s“^, which we subtract out. Further¬ 
more, the native CfA Digital Speedometer velocity sys¬ 
tem is s hifted by —0.14 km s“^ compared to the lAU 
system (jStefanik et al.l [19^ ^. The actual correction 
to get the TRES velocities onto the native CfA sys¬ 
tem is observed to be —0.75 km s“^, very close to 
—0.62 — 0.14 = —0.76 km s“^. We report TRES ve¬ 
locities on this native CfA system. 

We do not apply any shift to the WIYN RVs, and, 
as noted in Section IH we do not detect any significant 
zero-point shift between the WIYN RVs and the native 
CfA system. Therefore by construction, all RVs reported 
here, as well as the mean RV of the cluster, are on the 
native CfA system. 

3.1. Radial-Velocity Precision 

In order to detect binaries (and higher-order systems), 
we require multiple observations at multiple epochs with 
known RV precisions. Our sample contains observations 
of mainly narrow-lined stars within a modest magni¬ 
tude range, but taken at multiple different telescopes 
with different instruments. Therefore to facilitate our 
subsequent analysis of binarity (Section 15.IL we esti¬ 
mate single-measurement precision values for RVs de¬ 
rived from each respective telescope. 

The majority of the CfA measurements come from the 
MMT and the Tillinghast Reflector (both using the Dig- 
ital Speedometers) . We follow the same procedure as 
ICeller et al.l (j2008l) to empirically determine the typical 

^ Note that the offset of the CfA Digital Speedometer native 
velocity system reported in that paper has the wrong sign. 


single-measurement precision for observations made at 
these two telescopes and at the WIYN Observatory with 
Hydra. Specifically we fit a distribution of two degrees 
of freedom to the distribution of standard deviations of 
the first three RV measurements for each single-lined star 
in these samples, respectively. The results for all stars 
observed > 3 times at these telescopes are shown in Fig¬ 
ures [S] The dashed lines in each panel of Figure [3] show 
the best fitting x^ distribution functions, respectively, 
which yield typical single-measurement precision values 
of 0.8 km s“^ for RVs from the Tillinghast Reflector -|- 
DS, 0.7 km s“^ for RVs from the MMT, and 0.5 km s“^ 
for RV from WIYN. These single-measurement RV preci¬ 
sion values for the CfA telescopes are in good agreement 
with those found bv iMathieu et al.l (Il986l). a nd also agree 
with the value found bv iHole et al.l (120091) for observa¬ 
tions of NGC 6819 taken at the MMT and the Tilling- 
hast Reflector with the same instrumental setup. Like¬ 
wise, the single-measurement RV precision for observa¬ 
tions at WIYN found here is similar to that found for 
the narrow-lined stars in NGG 188, NGG 68 19 and M35 
observed with th is same setup (|Geller et al.ll200d 120101 : 
IHole et ahllMill . 

For these three telescopes (WIYN, MMT and Tilling¬ 
hast -I- DS), we also have sufficient data to characterize 
our RV precision as a function of V magnitude (which, 
in general, correlates with the signal-to-noise of the mea- 
surement for a give n int egration time). As al so found by 
iGeller et al.l (1200811 and iMathieu et al.l (|1986ll the single¬ 
measurement RV precision degrades towards fainter mag¬ 
nitudes. In Figure S] we show the single-measurement 
precision values resulting from x^ distribution fits to ob¬ 
servations from these three telescopes in bins of the stars’ 
V magnitudes. The lines show linear fits to these data. 
To avoid unrealistic precision values, we impose a preci¬ 
sion floor at the values in the bin containing the bright¬ 
est stars, specifically at 0.4 km s“^ for WIYN RVs and 
0.6 km s“^ for MMT and Tillinghast -|- DS RVs. We then 
use these fit lines to determine our single-measurement 
precision values for stars observed at these telescopes 
(given their respective V magnitudes) in later analyses. 

We also note that observations of the few rapidly ro- 
tating stars in o ur sample have poorer precision (see 
IGeller et al.ll20I0l for a detailed discussion of the effects 
of stellar rotation on the WIYN RV precision). For 
these few stars, we do not attempt to derive single¬ 
measurement RV precisions. 

We lack sufficient observations on the remaining tele¬ 
scopes to reliably utilize the x^ distribution fitting tech¬ 
nique. For the Palomar Hale 5m RVs, the typical stan¬ 
dard deviati on of the measurem ents of a given star is 
0.3 km s“^ (|Mathieu et al.l [19^ ; we use this value as 
our single-measurement RV precision for stars observed 
at Palomar. For observations from the Wyeth 1.5m, 
there are only 12 stars with at least three RV measure¬ 
ments, and a third of these stars are in binaries (with 
orbital solutions). The mean of the standard deviations 
of the first three RVs for the remaining eight stars is 
0.4 km s“^, and we take this as an estimate of the 
single-measurement RV precision for observations from 
the Wyeth. CORAVEL and TRES observations were 
primarily of binaries. For CORAVEL RVs, we estimate 
the RV precision based on the ( O — C ) residuals from 
the orbital fits in IMathieu et al.l (|I990f) . and find a typi- 
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Fig. 3.— Histograms of the RV standard deviations of the 
first three RV measurements for all stars observed > 3 times 
at the Tillinghast Reflector + DS (top), MMT (middle) and 
the WIYN (bottom) telescopes. Also shown in the dashed 
line in each panel is the best fitting function for each dis¬ 
tribution. The hts yield a single-measurement precision of 
0.8 km s“^ for Tillinghast + DS RVs, 0.7 km s“^ for MMT 
RVs and 0.5 km s"^ for WIYN RVs. See also Figure g] for 
an analysis of the RV precision as a function of the stars’ V 
magnitudes. 

cal precision of 0.5 km s“^. 

The correction of the TRES velocities for shifts in the 
instrumental zero-point from month to month and year 
to year i s now reliable at th e level of 0.02 km s“^ or better 
(e.g. see lQuinn et al.l[2014D . much better than the preci¬ 
sion achieved for most of our TRES spectra of M67 stars 
due to SNR limitations set by photon noise. Neverthe¬ 
less, most of our TRES observations of slowly-rotating 
single-lined stars in M67 do yield a velocity precison on 
the order of 0.1 (or perhaps 0.2) km s“^. We assume a 
single-measurement RV precision of 0.1 km s“^ for TRES 
observations of narrow-lined stars in the following anal¬ 



Fig. 4.— Single-measurement RV precision as a function 
of V magnitude for stars observed at the WIYN telescope 
(circles and solid line), MMT (triangles and dashed line) and 
the Tillinghast Reflector -|- DS (diamonds and dotted line). 
The points are placed at the centers of bins in V magnitude 
(with edges at V = 15.5, 14, 13, 12, 10), and show results 
from our fits to x^ distributions within these respective bins 
for each telescope. The lines show linear fits to these data, 
with a floor of 0.4 km s“^ for WIYN RVs and 0.6 km s~^ for 
both the MMT and Tillinghast -|- DS RVs. 

yses. 

TRES observations of rapidly-rotating stars have 
proven to be much more precise than the results for 
the same stars observed with the CfA Digital Speedome¬ 
ters. We attribute this to the higher SNR and better 
wavelength coverage provided by TRES. In particular, 
TRES has allowed us to derive reliable RVs for two of 
the rapidly rotating BSS in M67 that eluded success with 
our other facilities. 

The long tails in the distributions extending beyond 
the best htting curves shown in Figure [3] are populated 
by binary (and higher-order) systems. Moreover, because 
the standard deviation of the RVs from a binary should 
be significantly larger than the single-measurement RV 
precision, we use the results of this analysis to identify 
binaries in our sample. We discuss our criteria for iden¬ 
tifying binaries in Section [01 

4. THE COMBINED CFA AND WIYN DATA SET 

The Palomar velocities were adj usted to the nat i ve CfA 
velocity system as described in iMathieu et al.l (|1986ll . 
Here we check for a possible zero-point offset between 
the WIYN and CfA data by comparing the mean RVs for 
non-RV variable stars in both samples. Specifically we 
select only those stars that have three or more observa¬ 
tions in both samples, with P{x^) >0.01 and e/i< 3 (we 
discuss the derivation of these quantities in Section [5.ip . 
There are 95 such stars in our sample, and for each of 
these stars we calculate the mean RV from each sam¬ 
ple. The mean RVs for stars in the CfA sample that 
have observations at multiple telescopes are weighted by 
the respective precision values discussed above. We then 
calculate the star-by-star difference in the mean RVs, 
defined as RVcfA ~ RVwiyn, and find a mean difference 
of 0.008 km s“^, with a standard error of the mean of 
0.06 km s“^. As this difference is well below the RV 
precision of any telescope used here (and less than the 
magnitude of the offset of the native CfA system from 
the lAU system, as discussed in Section [3]), we conclude 
that there is no significant zero-point offset at this level 
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Fig. 5.— Histogram of the maximum time span of RV mea¬ 
surements for stars in our M67 sample that have at least 2 
observations. 


Figure [ 6 ] shows the percentage of stars with RV mea¬ 
surements in our stellar sample as functions of V magni¬ 
tude, {B — V) color and distance from the cluster center. 
There is no significant trend in our completeness within 
the primary sample with magnitude, color or radius. We 
show our completeness for stars with 15.5 < V < 16.5 in 
the left panel of Figure [5] within the gray region. There 
are 418 stars within this magnitude range that are within 
30 arcmin from the cluster center. We have at least one 
RV measurement for 295 (71%) of these stars, and >3 RV 
measurements for 182 (44%) of these stars. As shown in 
Figure [U the CfA sample also contains observations of 
stars located at >30 arcmin from the cluster center. We 
do not include these stars in our primary sample used 
in the subsequent analyses, but we include all observed 
stars in our RV summary table (Table IH presented in 
Section [5]). 

5. RESULTS 


between the WIYN and CfA RVs. 

We therefore proceed in combining the two data sets 
without modification, and use this combined data set for 
the following analyses. In total there are 8023 RVs of 
455 stars from the CfA and 5753 RVs of 1163 stars from 
WIYN for a total of 13776 RVs of 1278 stars in M67. 

In Figure [5] we plot a histogram of the maximum time 
span of RV measurements for the stars in our M67 sam¬ 
ple (with at least 2 observations). For some stars, our 
observations span >40 years, and about one third of the 
stars in our sample have observations spanning at least 
10 years. This long time baseline facilitates the detection 
of long-period binaries in our sample. Again, we discuss 
our binary detection method in Section [5.II 

4.1. Completeness of Spectroscopic Observations 

We define a “primary stellar sample” which extends 
from the brightest stars in the cluster to V < 15.5 and 
out to a 30 arcmin radius from the cluster center. In to¬ 
tal, there are 903 stars in our primary sample. We have 
at least one RV measurement for 902 of these stars, and 
>3 RV measurements for 899 of these stars. Moreover, 
there is only one star without any RV measurements; 
2038 (S2312)® is the second brightest star in our pri- 
mary sample, is brighter than our WIYN sample, and is 
a lSandei^ (Il977f) proper-motion non-member. The other 
three stars with < 3 RV measurements are 1016 (S1306), 
7033 (S1594), 7054 (F1295). 1016 is the brightest star 
in our sample and is a proper-motion non-member from 
three sources. We have two observations of 1016 that 
are outside of the cluster RV distribution. We have 
observed 7033 and 7054 multiple times, and they both 
appear to be rapid rotators. We have been unable to 
derive reliable RVs from the majority of these measure¬ 
ments. 7033 has one prop er-motion membership proba¬ 
bility ffrom lenders! HTTdl of 0%. 7054 has no proper- 
motion measurements. Our one tentative RV measure¬ 
ment of 7054 is outside of the cluster RV distribution, 
but is uncertain due to the rotation. In short, we have 
RVs for all likely cluster members within our complete 
sample. 

® See Section [ 5 ] for an explanation of our nomenclature for our 
IDs and cross-reference IDs. 


In the following sections we analyze the RV measure¬ 
ments of each star in our stellar sample. We first assess 
the RV variability of each star and then the membership 
probability. We use these quantities to classify stars with 
>3 observations as cluster members or non-members, 
and RV variables (e.g., binaries) or non-variables (e.g., 
single stars). We present the results of this analysis in Ta¬ 
ble IH where we provide the WOCS ID (/Dw)®, a cross- 
reference ID {IDxY, the J2000 right ascension {RA) and 
decli nation (Dec), the V magn itude and (B — V) color 
from I Montgomery et al.l (jl993l l. where available®, num¬ 
ber of RV measurements in the WIYN (Vw) and CfA 
(Nc) samples (counting RVs for primary stars only), the 
modified Julian date (JD - 2400000) of the first obser¬ 
vation (JDq) and last observation (JDf), the weighted 
mean RV (RV; for binaries with orbital solutions, we in¬ 
stead provide the center-of-mass, 7 -RV) and the weighted 
standard error of the mean RV (RI4; for binaries with or¬ 
bital solutions we instead provide the error on the q-RV), 
the c ombined RV precision (i, defined as in iHole et ^ 
1200911 the RV membership probability (Prv, defined in 
Section 15.2L the prop e r-mot ion me mbership probabili- 
ties from Yadav et al.l (l2008l PpMy), IZhao et Si ( 1993 . 
PpMz), iGirard et al.l (|l989l and ISandersI ( 1977 . 

PpMs), where available, the e/i and P{x^) values (de¬ 
fined in Section EH), the membership classification (see 
Section [5.3p and finally a comment field. 

We mark particularly notable stars within our stel- 


® We follow the same method as IHole et aLI l|2Q09ll to define 
WOCS ID based on the given star’s V magnitude and distance 
from the cluster center. 

^ If available the cross-reference ID is tak en fromISanderd (119771) 
and denoted by the pre fix “S”. If there is n olSanderd (|1977l ) source 
match, we provide the IMontgomerv et al.l l|1993D ID. if available, 
denoted by the prefix “M”. If bo t h of t hose studies lack a source 
match, we provide the lFan et al.l 19961) IDs, if available, denoted 
by the prefix “F”. There are 13 sources in our table that do not have 
matches in these three references. As with all sources, we provide 
their RA and Dec positions for matching to other catalogs. 

^ For the few stars in our sample that are not in the 
IMontgomerv et al.l l| 19931) survey and for which we can find no 
other BV photometry from the literature, we derive V magni¬ 
tudes from the 2MASS JHK photometry using a similar rela¬ 
tionship to IGirard et al.l (|2QQ4I ). For these stars, we do not have 
[B — V) colors. The stars wi t h IDs that have a prefix of “T” 
are taken from IMathi eu et al.l (|1986l) who use photometry from 
IMurrav &: Clements! (1196^ 
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Fig. 6. — Completeness histograms as a function of apparent V magnitude (left), (B — V) color (middle) and radius from the 
cluster center (right) for stars observed in M67. In the V completeness panel we show all stars observed in this survey. The 
gray region shows stars that are fainter than our primary sample selection of L < 15.5. (Note that there are only two stars 
brighter than V = 8.5 in our sample, both of which are proper-motion non-members and neither of which have >3 observations, 
as discussed in the main text.) In the (B — V) and radius panels we show our completeness for stars in this primary sample 
only. In all plots the percentage of stars observed >1 time is plotted with the solid line, while the percentage of stars observed 
>3 times is plotted with the dotted line. 


lar sample in this c omment field of Ta ble ID X-ray 
sources identified by iBelloni et al.l (|1998D with ROSAT 
are labeled with “X” followed by the source number 
giv en in their paper. Likewise X -ray sources identified 
bv IVandenBerg fc StetsonI (I2004D with Chandra are la¬ 
beled with “CX” followed by the sourc e number given 


I99II: iStassun et al.l 120021 

Ivan den Berg et al.l 

2002 

Sandouist & Shetronel l2003allbl: iSandouist et al.l 

2001 

Oian et al.1120061: iBruntt et al.l 

20071: iStello et al.l 

20071; 

Pribulla et al.1 20081: [Yakut et al.ll2009D are labeled with 


“PV”, or “PV?” if the authors identify the photometric 
variability as uncertain (e.g., possible flare event detec¬ 
tions, possible nearby source contamination, etc.). We 
also label the W UMa’s found in M67, and provide the 
GCVS names for photometric variables, where available 
(e.g., AH Cnc, ES Cnc, etc.). Stars that are rotating 
significantly more rapidly than our instrumental resolu¬ 
tion of ^10 km s“^ are labeled as “RR”. Sources that 
are detected as triple systems are labeled with “triple”. 
Additionally, we label blue stragglers with “BSS”, “yel¬ 
low giants” with “YG”, and the two sub-subgiants with 
“SSG”. We briefly discuss a few of these notable stellar 
populations in Section [Hill 


5.1. Radial-Velocity Variability 

We identify binaries in our M67 data using the eji 
statistic, which is the ratio of the standard deviations 
(e) to the expected precision (i) of the RVs for a given 
star. Members of binaries with orbital periods short 
enough to show significant RV variation in our data will 
have higher standard deviations than expected for sin¬ 
gle stars. Therefore high eji values indicate binaries or 
higher-order systems. As a large number of stars in our 
sample have RV measurements from multiple telescopes, 
and therefor e multiple precision values, w e use the for¬ 
malism from iBevington fc R.obinsonl (Il992f l to derive the 
e and i values for da ta with multiple precision values, as 
in iHole et al.l (l200^ . 

We derive aneji value for all single -lined stars with >3 
RV measurements. Previous work bv iGeller et al.l (|2008[1 
shows that stars with e/i>4 can be securely identified as 
binaries (or higher-order systems), and we initially fol¬ 


lowed this same cutoff for identifying binaries here. In 
practice, however, we have derived RV orbital solutions 
for all proper-motion members in our primary sample 
with 3 < e/i < 4. (There are only three proper-motion 
non-members in our primary sample with 3 < e/z < 4; 
each of these have mean RVs outside of the cluster distri¬ 
bution.) Therefore here, we identify binaries in our data 
as having e/z>3. 

The uncertainties for RVs in double-lined (SB2) bina¬ 
ries are not well established, and therefore we do not de¬ 
rive eji values for these stars. Instead we identify SB2s 
as binaries directly by inspection of the spectra and the 
morphology of the peaks of the cross-correlation func¬ 
tions. (We derive the uncertainties on the mean RVs, 
i?I4, for SB2s without orbital solutions using the mea¬ 
surement precisions that we derive above for single-lined 
stars; we suspect that the RV uncertainties for such SB2s 
quoted here are in fact lower limits.) 

In addition to the eji statistic, we also provide the 
P{x^) value. This statistic tests the hypothesis that a 
given star’s distribution of RVs is consistent with a con¬ 
stant value at the mean RV. Members of binaries with 
short enough orbital periods to show significant RV vari¬ 
ations in our data will be inconsistent with a constant RV 
and therefore will have small P{y^) values. To derive the 
P{x^) value, we first calculate the standard statistic, 
using the weighted mean RV and the respective precision 
of each RV measurement (discussed in Section [XT]) . For 
each observed star, we use the number of RV measure¬ 
ments minus one for the degrees of freedom. Then P{x^) 
is the corresponding probability for obtaining a value of 
X^ greater than or equal to the observed value with the 
given degrees of freedom. Again, we only derive P(x^) 
values for single-lined stars with >3 measurements. 

In the following analysis we will identify binaries by 
having e/i> 3, or having a binary orbital solution. (The 
additional P{x^) statistic is provided for interested read¬ 
ers who prefer to use these values to perform their own 
selection of binaries.) All stars with e/z < 3 for which we 
have not derived binary orbital solutions, are labeled as 
“single”. However, some of these stars are undoubtedly 
in long-period binaries, currently beyond our detection 
limit. 
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TABLE 2 

Gaussian Fit Parameters For Cluster and Field 
Radial-Velocity Distributions 


Parameter 

Cluster 

Field 

Ampl. (Number) 

147.8 ± 0.7 

2.38 ± 0.10 

RV (km s-l) 

33.615 ± 0.005 

23.4 ± 2.2 

a (km s“^) 

0.854 ± 0.005 

46.7 ± 2.2 


Finally, we also identify five likely triple stars in our 
sample, and label them as such in TablelH In most cases 
the systems are double lined, where the RVs for one star 
vary on a much shorter time scale than those of the other 
star. For a few, we see hints of tertiary velocities at low 
signal-to-noise, but we have yet to analyze these spec¬ 
tra in detail to derive all three velocities simultaneously. 
Additional triple stars may be detectable in our sample, 
for instance within the RV residuals from binary orbital 
solutions, but we save this analysis for a future paper. 


5.2. Membership 

In Figure [3 we show the distribution of weighted mean 
RVs for all stars in our sample with >3 RV measure¬ 
ments and e/z < 3, and y-RVs for binaries with orbital 
solutions. For stars with RVs from multiple telescopes, 
we use the respective RV precision values to calculate the 
weighted means, and use those results here. The cluster 
RV distribution is clearly distinguished from that of the 
held as the narrow distribution peaked at a mean RV of 
-1-33.6 km s“^. 

In order to derive membership probabilities, we hrst 
ht simultaneous one-dimensional Gaussian functions to 
the cluster and held RV distributions, Fc{v) and Ff{v), 
respectively. (We exclude from this ht binaries without 
orbital solutions, as their y-RVs are unknown.) The re¬ 
sulting combined ht to the cluster and held distributions 
are shown in the blue line in the left panel of Figure [3 
and the ht parameters are given in Table [2] (ht to a his¬ 
togram with a bin size of 0.5 km s“^). We then use the 
following equation. 


Prv{v) 


Fcjv) 

Ff{v)+Fc{v) 


( 1 ) 


(|Vasilevskis et al.lll958D to calculate the RV membership 
probability Pry(v) for a given star in our sample. 

We only compute membership probabilities for stars 
with > 3 observations. For non-RV-variable stars, we use 
the weighted mean RV in the calculations. For binaries 
with orbital solutions, we use the y-RV. For RV-variable 
stars without orbital solutions, we cannot calculate a re¬ 
liable RV membership probability, as the y-RV is un¬ 
known. For these stars we instead provide a preliminary 
membership classihcation, described in Section [573] 

We plot the resulting distribution of RV membership 
probabilities in the right panel of Figure [3 Cluster and 
held stars are cleanly separated, and we choose a cutoff 
of Trv> 50% to dehne our M67 cluster member sample. 
Using our sample of single cluster members with e/z < 3 
and binary cluster members with orbital solutions, we 
hnd a mean cluster velocity of -1-33.64 km s“^ (with an 
internal precision of ±0.03 km s“^), in good agreement 
with previous RV surveys (e.g., see lYadav et al.l [20081 
and references therein). 


As stated in Sections [3] and |4l the RVs in this paper 
are all on the native CfA system, which is found to be 
shifted by —0.14 km s“^ compared to the lAU system. 
Therefore, the mean cluster velocity quoted here may 
also have this offset compared to the lAU system (i.e., 
±0.14 km s“^ may need to be added to our velocities to 
get onto the lAU system). 

The ratio of the areas under the Gaussian hts to the 
cluster and held distributions provides an estimate of the 
held star contamination. At a membership probability of 
50%, we expect a 3.5% contamination from held stars in 
our RV cluster member sample (i.e., ^ 20 stars) . 


5.2.1. Comparison of Radial-Velocity and Proper-Motion 
Membership Probabilities 

Here we compare with four pu blished proper - 
motion menibership catalogs fo r M67 : ISar 
iGirard et al.l (Il989li . IZhao et all (|1993ll and lYadav et al.l 
( 200811 . For this comparison we use only non-RV variable 
stars and binaries with orbital solutions in our primary 
sample (thereby excluding binaries whose 7 -RVs are un¬ 
known) so as to ensure secure RV membership probabil¬ 
ities. There are 836 stars in our primary sample that 
meet these criteria. 

We hnd excellent agreement with all four proper- 

moti on sourc es. _ There are 456 stars in this sample 

with ISander j (|1977ll proper-motion membership proba¬ 
bilities >50, of which 430 (94%) also have RV member- 
ship probabi li ties > 50%. 410 stars in this sample have 
IGirard et al.l (Il989| l proper-motion membership proba¬ 
bilities >50%, of which 396 (97%) also have RV member- 
ship probabilities >50%. 164 stars in this sample have 
IZhao et al.l (|1993ll proper-motion membership probabil¬ 
ities >50%, of which 149 (91%) also have RV member¬ 
ship pro babilities >50%. F inally, 330 stars in this sam¬ 
ple have lYadav et all (1200811 proper-motion membership 
probabilities >50%, of which 200 (91%) also have RV 
membership probabilities >50%. If we examine stars in 
this sample that have proper-motion membership prob¬ 
abilities of >50% from all four references, we find that 
95/100 (95%) also have RV memberships >50%. 

For stars with both RV and proper-motion membership 
probabilities, we combine the results to further refine our 
cluster member sample. We first take all stars that have 
RV membership probabilities >50% as candidate clus¬ 
ter members. We then examine the available proper- 
motion membership probabilities for each star to remove 
non-members from this sample. For each proper-motion 
reference we assign the same cutoff of >50% member¬ 
ship probability to denote cluster members. If the avail¬ 
able proper-motion membership values for a given star 
all indicate that the star is a non-member, we remove 
it from our cluster member sample. For stars with 
proper-motion membership probabilities from multiple 
references that disagree on membership status, we allow 
each proper-motion reference one “vote” on membership 
(either member or non-member) and take the majority 
vote to determine the proper-motion membership sta¬ 
tus. If this procedure results in a tie, we instead use the 
result from the highest precision proper-motion member¬ 
ship survey for t he given star to determine membership 
(in the order of IGirard et al.l . 9891 lYadav et al.l l2008l 
Zhao et aI]ll993L thenISanders I197 a with IGirard et al.l 


19891 at the highest precision). 
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Fig. 7.— RV distribution (left) and the distribution of RV memberships (right) for stars in our M67 stellar sample. To derive 
the RV distribution we use the weighted mean RVs for stars with >3 observations with e/i< 3, and 7 -RVs for binaries with 
orbital solutions. (Here we exclude detected binaries without orbital solutions, because their 7 -RVs are unknown.) We use 
bins of 0.5 km s“^. The red line shows the field component to our simultaneous two-Gaussian fit to the observations, and the 
blue line shows the combined fit of the cluster and field distributions, which we use to derive the membership probabilities for 
observed stars in M67. On the right, we show a histogram of RV memberships (rounded to the nearest percent) of the same 
stars used in constructing the RV distribution. RV members are defined as having Prv> 50%. 


In practice, only 24 stars in our primary sample have 
proper-motion membership probabilities that result in a 
“tie”, all but two of which have only two proper-motion 
membership values. Also for reference, there are 35 stars 
in our primary sample with proper-motion memberships 
from all four studies where one disagrees wit h the other 
three. In 23 of these ca ses. IZhao et al.l (Il993ll is the out¬ 
lier. In 8 of these c ases. lYadav et all (l2008l) is the outlier. 
In 4 of these cases. [Girard et ahT i I989D Is the outlier. 

Of the 515 stars in our primary sample that have se¬ 
cure RV memberships >50% (i.e., SM and BM stars, see 
Section 15.311 and proper-motion measurements, 455 pass 
this proper-motion membership test, and are therefore 
deemed bona fide cluster members. This ^12% contami¬ 
nation of our RV member sample is about 3.5 times larger 
than we estimate above based solely on the RV distribu¬ 
tions of the cluster and field stars. We note that 75% 
of these RV members that we determine to be proper- 
motion non-members from the algorithm above reside 
outside of 15 arcmin from the cluster center. Also about 
23% of these RV members that are proper-motion non¬ 
members have at least one proper-motion membership 
value (from one of the references above) greater than 
50%. However, to minimize possible field star contamina¬ 
tion, we remove all of these proper-motion non-members 
from our cluster member sample. 

5.3. Membership Classification of Radial- Veloeity 
Variable Stars 

For stars that show no significant RV variability and 
binaries with orbital solutions, we can calculate precise 
RV membership probabilities that allow us to separate 
cluster members from field stars, as described above. 
However, for binary stars that do not have RV orbital 
solutions, we cannot calculate reliable RV membership 
probabilities because we do no t know their 7 -RVs . Thus 
we f ollow a sirn i lar m ethod to iGeller et al.l (j2008L 120101) 
and iHole et all (|2009D in order to provide a qualitative 
classification of the membership and variability for each 
observed star. 

Our classification scheme is defined in Table O where 
we identify our eight qualitative membership classes, the 


selection criteria for each class, and also the number of 
stars that reside in each class. The selection criteria de¬ 
pend on the number of RVs, the eji value, the RV mem¬ 
bership probability (Prv, determined either using the 
mean RV, RV, or the 7 -RV, for binaries with orbital so¬ 
lutions), and the proper-motion membership “vote” de¬ 
scribed in Section [5. 2.1 1 For the proper-motion member¬ 
ship vote, we use “PpM=]VI” to indicate proper-motion 
members, and “PpM=NM” to indicate proper-motion 
non-members. If a star does not have a proper-motion 
measurement, we use only the RV criteria to classify the 
star. 

In short the single members/non-members (SM/SN) or 
binary members/non-members (BM/BN) are stars with 
secure membership status; these are the only stars for 
which we can provide reliable RV membership probabil¬ 
ities. For RV variable stars with >3 RV measurements 
that do not have a binary orbital solution, we divide our 
classification into three groups, including binary likely 
members (BLM), binaries with unknown RV membership 
(BU) and binary likely non-members (BLN). We antic¬ 
ipate that eventually orbital solutions derived for BLM 
binaries will place them in the BM category, while those 
for the BLN binaries will place them in the BN cate¬ 
gory (since many of these sources are proper-motion non¬ 
members, and for those that are not, it is unlikely that an 
orbital solutions will place their 7 -RVs within the cluster 
distribution). Binaries with unknown RV membership 
(BU’s) are proper-motion members. Therefore here we 
assume that these are indeed cluster members (unlike in 
other WOGS papers, where proper-motion memberships 
were unavailable). Stars with <3 RV measurements have 
unknown RV membership (U), as these stars do not meet 
our minimum criterion for deriving RV memberships or 
eji values. In the following analyses we restrict our clus¬ 
ter member sample to only include the 562 stars classified 
as either SM, BM, BLM or BU. 

Finally, as mentioned above, stars that show broad¬ 
ened spectral features (e.g., due to rapid rotation) do 
not have secure single-measurement RV precision values, 
and therefore in most cases we cannot confidently clas¬ 
sify such stars as binaries or singles. We provide our best 
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TABLE 3 

Description of Stars or Star Systems Within Each Membership Class 


Class 

Description 

Criteria 

Number 

SM 

Single Member 

>3 RVs, e/i < 3, ^ 50% AND PpM— 1'^ 

420 

SN 

Single Non-member 

>3 RVs , e/i < 3, < 50% OR iV’M— NM 

498 

BM 

Binary Member 

binary orbit, Prv( 7 -RV) > 50% AND -PpM= M 

108 

BN 

Binary Non-member 

binary orbit, fRv( 7 -RV) < 50% OR PpM= NM 

17 

BLM 

Binary Likely Member 

>3 RVs , e/z > 3, Prv(RV) > 50% AND PpM= M 

23 

BU 

Binary with Unknown RV 
Membership 

>3 RVs, e/i > 3, Prv(RV) < 50% AND the range in RV measurements includes 
the cluster mean RV AND PpM= M 

11 

BLN 

Binary Likely Non-Member 

>3 RVs, e/i > 3, Prv(RV) < 50% AND all RV measurements are either at higher 
or lower RV than the cluster distribution, OR PpM= NM 

58 

U 

Unknown RV Membership 

< 3 RVs 

143 


assessment of the binarity of these sources in Table [Hand 
indicate our uncertainty in their class with parentheses, 
e.g. (BL)M, (S)N, etc. 

6 . DISCUSSION 

In the following section, we use our confirmed clus¬ 
ter members to investigate the CMD (Section [RT]), iden¬ 
tify and discuss a few notable stellar populations (Sec¬ 
tion |621), analyze the spatial distribution of the single, 
binary, giant and BSS cluster populations fSection 16.3|) 
and derive the velocity dispersion of the solar-type stars 
in our cluster member sample ISection 16.41) . 

6 .1. Color-Magnitude Diagram 

In Figure [5] we plot the CMD of all stars in our stel¬ 
lar sample with V < 15.5 (and available {B — V) colors; 
left) and only the confirmed cluster members within the 
same magnitude range (right). Without removing non¬ 
members, the main-sequence of the cluster is visible, but 
the BSS and giant populations cannot be distinguished 
from the field. Applying both our proper-motion and RV 
membership criteria reveals a rich cluster containing well 
populated main-sequence, subgiant and giant branches 
as well as a large population of BSS (blue points), four 
yellow giants (red points), and two sub-subgiants (green 
points). Specihcally, we plot here stars that we clas¬ 
sify as SM, BM, BLM and BU, and we also include here 
stars in the U category as well as rapid rotators that are 
proper-motion members. 

For comparison, we also plot a 4 Gyr isochrone and a 
zero-age main-sequence isochrone (solid lines), as well as 
an equal -mass binary l i ne (da shed line). The isochrones 
are from iMarigo et all (|2008ll . and use solar metallicity, 
(m — M)v = 9.6 and E{B — V) = 0.01, consistent with 
recent results derived in the literature (see Section [T]). 
We include the isochrones simply to help guide the eye; 
they are not meant as a ht to the observed data. 

Binaries with orbital solutions are circled, and binaries 
without orbital solutions are marked with diamonds. As 
is seen clea rly here, a nd als o noted by iLatham fc Miloii^ 
(|l996ll and iLathaml (|2007ll . the “exotic” stars (i.e., the 
BSS, yellow giants, and sub-subgiants) have a remark¬ 
ably high binary frequency. In total we identify 14 BSS, 
four yellow giants (one of which is outside of a 30 arcmin 
radius from the cluster center) and two sub-subgiants in 
our stellar sample. At least 16/20 (80% ± 20%) of these 
exotic stars are RV variables (and others may also be 
binaries, e.g. with long-period orbits that are currently 


outside of our detection limit). In comparison, 122/538 
(22.7% ± 2.1%) of the “normal” stars, located in more 
typical positions in the CMD, show RV variability indica¬ 
tive of binarity. Thus the exotic stars have a signihcantly 
higher frequency of binaries than the normal stars. This 
result is similar to that found in the old (7 Gyr ) open 
cluster NGC 188, where iMathieii fc GelleJ (|2009l l2015ll 
find that 80% of the NGC 188 BSS have binary com¬ 
panions (roughly three times the binary frequency of the 
main-sequence stars in NGC 188). 

We also note the very tight red giant seq uence, first dis¬ 
cussed in detail bv Danes fc SmitlJ (|1984ll . Such a tight 
red giant sequence is expected for a single coeval pop- 
ulation, but inter estingly, the giants of NGC 188 (e.g. 
iGeller et al.l 1200^ and the in termediate-age ( 2.5 Gyr) 
open cluster NGC 6819 (e.g. iHole et all l2009ll show a 
much larger scatter than the giants in M67. The origin 
of the scatter in these other open clusters is unknown. 

6 .2. Stars of Note 

Within our cluster member sample, there are a num¬ 
ber of intriguing stellar populations, from those that lie 
far from the predicted locus of single stars from stel¬ 
lar evolution theory (including the well known BSS, Sec- 
tion l6.2.11 and anomalous giants, Section lB.2.21) . to those 
that have physical characteristics very similar to those of 
our Sun ISection I6.2.3p . and even three exoplanet host 
stars ISection 16.2.41) . In the following we briefly identify 
and discuss these stars of note. 

6.2.1. Blue Stragglers 

M67 is home to one of the most well-studied BSS 


popu lat ions (see, e.g ., lEggenl Il981t IPeterson et al.l 


19841: 

Mathieu et al. Il986l iManteiga et al.l 19891 


Gilliland & BrownI 119921: iLeonard & Linnelll 

19921: 

Vlilonel 119921: iMilone & LathamI 11994 iLeonardI 


I/andsman et al.l 11997, 119981: iDeng et al.l 

1999t 

Shetrone & Sandauisd 20001 iHurlev et al. 


2 ^ 

van den BerEf et al.l l2001t iSandauist et al. 

2003t 

Chen & HanI 

2004 iHurlev et al. 120051 Sandauistl 

2005 

Zhang et al.l 

20051: Andronov et al.l 120061: Tian et al.l 

20061 iBruntt et al.ll2007t ILathaml 120071 iLiu et al.l 120081: 

Pribulla et al.l 20081: iLu et ^l.|2010|l. In total, we identifv 


25 candidate BSS that are discussed in the literature: 


1006 (S1066), 1007 (S1284), 1010 (S977), 1017 (S1466), 
1020 (S751), 1025 (S1195), 1026 (S1434), 2007 (S984), 
2008 (S1072), 2009 (S1082), 2011 (S968), 2013 (S1267), 
2015 (S792), 2068 (S277), 3005 (S1263), 3009 (S1273), 
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Fig. 8 .— Color-magnitude diagrams for all stars in our M67 stellar sample with V < 15.5 (and available {B — V) colors; left) 
and only cluster members within the same magnitude range (right). We take here as cluster members all stars that reside in 
the SM, BM, BLM, and BU classes, as well as stars in the U class and rapid rotators that are proper-motion members (see 
Section [5.31) . Binary members with orbital solutions (BMs) are circled, and velocity variables without orbital solutions that are 
likely cluster members (BLMs and BUs) are shown in dia monds. For compar ison, in the right panel we also plot a zero-age 
main-sequence isochrone and a 4 Gyr isochrone (solid lines iMarieo et ^120081 ') using (m — M)v = 9.6, E{B — V) = 0.01 and 
solar metallicity. We also plot the equal-mass binary locus (dashed line) obtained by shifting the 4 Gyr isochrone by -0.75 mag. 


3010 (S975), 3013 (S752), 4003 (S1036), 4006 (S1280), 
5005 (S997), 5071 (S145), 6038 (S2226), 8006 (S2204), 
9005 (S1005). 

We find 1017 and 5071 to be cluster non-members 
based on kinematic membership information, and there¬ 
fore we remove these from the list of BSS. 1017 is a non¬ 
member by both proper motions and RVs. 5071 has 
a prope r -motio n membership probability of 94% from 
LSandersI (|1977[ 1. but the star appears to be single with 
a mean RV or 39.79±0.12 km s“^, which results in a 0% 
RV membership probability. 

By examination of the CMD shown in Figure [U we 
remove an additional nine stars from the BSS sample. 
1020, 4003, 8006 and 9005 lie close to the bl ue hook of the 
cluster (given the iMontgomerv et al.lll99^ BV photome¬ 
try), and therefore we conservatively exclude these stars 
from our BSS sample. 2007, 2015, 2068 and 3009 are 
above the turnoff, but reside in a region expected to be 
populated by binaries containing normal main-sequence 
turnoff stars. We only detect a binary companion to 
2068. However, the remaining three stars may have long- 
period companions, currently beyond our detection limit, 
with high enough masses to explain their location on the 
CMD. We therefore do not include these stars in our BSS 
sample. 2008 is the reddest BSS candidate in this sam¬ 
ple, residing ^1.5 magnitudes directly brighter than the 
main-sequence turnoff. The BSS status of 2008 is some¬ 
what ambiguous, as it may currently be evolving toward 
the giant branch. We choose to exclude 2008 from our 
BSS sample and will inclu de it in our “anomalous gi ant” 
sample discussed below. (jMathieu fc Lathirinl Il986l also 
exclude 2008 from their BSS sample based on similar ar¬ 
guments.) 


Thus in total we find M67 to have 14 BSS in our pri¬ 
mary stellar sample. Remarkably, 11 of these 14 BSS 
show significant RV variations indicative of binary com¬ 
panions. Seven of these BSS have secure orbital solu¬ 
tions, and an additional two rapidly rotating BSS (1026, 
4006) have preliminary orbital solutions. The detected 
binary frequency amongst the M67 BSS is 79% ± 24%. 
This high binary frequency for the M67 BSS is similar 
to that found in the old (7 Gyr) open clu ster NGC 188 
(iMathieu fc Gelle'ill2009D . and in the field (|Carnev et al.l 
I 2 OOII) . We will discuss these BSS in detail, including 
their binary properties, in a subsequent paper. 

6.2.2. Anomalous Giant Stars 

There are six stars with (B — V) colors consistent 
with the subgiant or giant branches but have magnitudes 
and/or chemical abundances that set the apart from the 
more typical giants along the isochrone shown in Fig¬ 
ure m We discuss these anomalous giant stars below. 

Sub-subgiants: 15028 (S1113) and 13008 (S1063) 
are fainter than the subgiant branch but redder than 
the main-sequence. These sub-subgiants are both mem¬ 
bers of binary systems and are also X-ray sources. 
Their origins are unk nown, and we refer the reader to 
IMathieu et al.l (120031) for a very detailed discussion on 
the available observations for these stars. _ 

Lith ium-rich Subgiant: iGanto Martins et al.l 

(|2006t) find the subgiant 6008 (S1242) to have an 
anomalously high Lithium abundance as compared 
to normal main-sequence turnoff stars in the cluster, 
although it falls within the normal subgiant branch on 
the GMD. We confirm the membership of 6008, and we 
also confirm that 6008 is in a binary. 
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Yellow Giants: 1015 (S1237), 1112, 2002 (S1040) and 
2008 (S1072) are all brighter than the giant branch and 
have often been referred to as “yellow giants”. 1112 is 
55.6 arcinin from the cluster center, and is therefore not 
included in our primary stellar sample. This star shows 
no evidence for a binary companion. The other three yel¬ 
low giants a re all members of binaries. 2002 was studied 
in detail by iLandsman et al.l (|1997L Il998ll who find the 
secondary to be a low-mass He white dwarf, suggesting 
that the system went through an episode of mass trans¬ 
fer while the donor was on the giant branch. Therefore 
2002 may have been a BSS in the recent past, and is 
now in the process of evolving towards the giant branch. 
2008 may also have previously bee n a BSS that is now 
evolv ing towards the giant branch (jMathieu &: LathamI 
1198611 . A similar evolutionary scenario may explain the 
anomalous CMD location of 1015. 

6.2.3. Solar Twins 

M67 is of a similar age and chemical composi¬ 
tion to the Sun, and is therefore an ideal tar¬ 
get for investigating solar analogs (iPasauini et ahl 
l2008t iReiners fc GiamDaDall2009t iGiampana et al.1120061 : 
ICastro et al.l 120111 : lOnehag et al.l l201lh . Stars 7003 
(S1041), 10018 (S1462), 11018 (S1095), 12012 (S996), 
13021 (S945), 14014 (S779), 16011 (S770), 16023 
(S2211), 17026 (S1335), and 18013 (S785) were found 
to have effective temperatures consistent with the Sun, 
and are the closest analogs to the Sun in M67. We con¬ 
firm that these 10 stars are cluster members, and all but 
one appear to be single (with 10018 at e/z =4.19). 16023 
was fo und to be a photometric variable bv IStassun et al.l 
(j2002tl , although variability was only detected above the 
3(7 level in the B band (and not in V or J). We find 
16023 to have e/i=0.93. 

6.2.4. Exoplanet Host Stars 

iBriicalassi et al.l (|2014l l identify three stars hosting 
roughly Jupiter-mass exoplanets in M67: 1045 (S364), 
13014 (S802) and 16011 (S770). Our measurements are 
not sensitive enough to detect such low-mass compan¬ 
ions, but we confirm that none appear to have stellar- 
mass companions and that all three are cluster mem¬ 
bers. ^We_also_j^e_Hmt 1045 is an X-ray source (X19 
from iBelloni et al.lIl998| l . 

6.3. Radial Distribution of Cluster Members 

The spatial dist ribution of stars in M67 has been stud¬ 
ied in detail (e.g. iMathieu fc LathamI 1198'^ IZhao et al.1 
1996t ISaraiedini et al . Il999t iBonatto fc Bical l2003l 
Davenport fc Sandauistl[2?)10l l. M67 is mass segregated, 
as is expected for a 4 Gy r cluster with a half-mas s re¬ 
laxation time of 100 Myr (jMathieu fc Lathamlll986ll . In 
light of our new RV memberships and identification of 
binaries, we briefly re-investigate the spatial distribution 
of cluster member populations using our primary sam¬ 
ple. In Figure[9]we show the projected radial distribution 
for the single main-sequence stars, binary main-sequence 
stars, single giants and BSS. Singles are stars classified 
as SM, and binaries are classified as either BM, BLM or 
BU. 

A Kolmogorov-Smirnov (K-S) test shows that the 
main-sequence binaries are centrally concentrated with 
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Fig. 9. — Cumulative projected radial distribution of cluster 
members in M67. We divide our primary sample into single 
and binary main-sequence (MS) stars, single giants and blue 
stragglers. The giants are identified as having {B — V) > 0.9, 
and we exclude the anomalous giants discussed in Section[6]2] 
The blue straggler sample is identified as discussed in Sec¬ 
tion lO and shown in Figure jS] Both the binaries and blue 
stragglers are significantly centrally concentrated with respect 
to the single main-sequence stars. 

respect to the main-sequence single stars at the 98.7% 
confidence level. Comparing all single stars to all bina¬ 
ries shows that the binaries are more centrally concen¬ 
trated at the 99.8% confidence level. Becanse the total 
mass of a given binary is generally more massive than 
that of a single star in our primary sample, this results 
c onhrms that M67 is mass se gregated. 

IMathieu fc LathamI (|I986tl find the BSS to be cen¬ 
trally concentrated with respect to the single stars near 
the cluster turno ff. Our BSS sample is sli ghtly differ¬ 
ent than that of IMathieu fc LathamI (|1986f ). Nonethe¬ 
less, the result is the same. A K-S test comparing the 
BSS to the single main-sequence stars show that the BSS 
are centrally concentrated at the 99.7% confidence level. 
This result suggests that the BSS are more massive than 
the main-sequence stars in the cluster, as was noted by 
IMathieu fc Lathan] (II986D . 

In many globular clusters and also the old open cluster 
NGC 188, BSS show a bimodal projected radial distri¬ 
bution (e.g.lFerraro et al.fl997Ll2012l : lMapelli et al.l[2006l : 
iGeller et al.ll2008l l. We do not observe evidence for a bi¬ 

modal radial distribution in our M67 BSS sample. How¬ 
ever, our sample only extends to 6 or 7 core radii. In 
NGC 188, which is of a similar dynamical age to M67, 
the halo BSS population begins at roughly 6 core radii 
and extends to about 13 core radii (at least). Thus if M67 
and NGC 188 have similar BSS radial distributions, we 
would only expect to see the inner population of M67 in 
our current sample, and would require a survey extend¬ 
ing to roughly twice the current radial extent to search 
for a bimodal structure. 

Interestingly, the giants appear to follow a similar spa¬ 
tial distribution as the binaries, despite having very sim¬ 
ilar masses to the upper main-sequence stars included in 
our primary sample. However, given the relatively small 
sample size of the single giants, a K-S test comparing the 
spatial distributions of the single giants and single main- 
sequence stars returns a distinction at only the 95.9% 
confidence level. For comparison, the NGC 188 giants 
do not follow the distribution of the binaries, and instead 
follow closely to the single cluster members (see Figure 9 
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in iGeller et alJ[2008f) . both when considering the entire 
spatial extent of the RV survey and when limiting to a 
similar number or core radii as our M67 sample. This 
difference is intriguing, but the small sample sizes of the 
giants in both clusters make such comparisons uncertain. 

6.4. Cluster Radial- Velocity Dispersion and Virial Mass 

We follow the method of iGeller et al.l (|2010t) to calcu¬ 
late the RV dispersion of the M67 single main-sequence, 
subgiant and giant members in our primary sample. (We 
exclude the binaries and BSS because these have signif¬ 
icantly different spatial distributions, likely due to their 
higher masses.) This method assumes that the observed 
RV dispersion is composed of two components, one which 
we will call the “combined RV dispersion” and one for 
the observational error, and also assumes that the distri¬ 
butions of RVs and errors are Gaussian. The observed 
dispersion is directly measured, and by subtracting off 
the component from observational error, we recover the 
combined RV dispersion. The combined RV dispersion 
itself contains contributions from both the true RV dis¬ 
persion of the cluster and undetected binaries, which also 
artificially inflates the observed RV dispersion. We aim 
to recover the true RV dispersion from our data. 

First we calculate the observed RV dispersion, which 
is simply the standard deviation of the observed mean 
RVs for each star about the cluster mean RV (calcu¬ 
lated in Section 15.21) . In our calculation, we use the 
weighted mean RVs for each star in the given sample 
(see Table S]) , which utilize the single-measurement pre¬ 
cision values from the given telescope for each individ¬ 
ual RV value (see Section [XT]) . For the single main- 
sequence, subgiant and giant stars in our primary sam¬ 
ple, the resulting observed velocity dispersion is 0.86 ± 
0.11 km s“^. After correcti ng for the observational er- 


ror (following the method of 

McNamara & Sanderslll977l 

and iMcNamara & Sekignchi 

1198611. we find a combined 


RV dispersion of 0.80 ± 0.04 km s“^. 

We then follow a similar method to that of iGeller et al.l 
(|2010ll to correct this combined RV dispersion for the 
contribution from undetected binaries. Briefly, we run a 
Monte Carlo analysis to create many realizations of our 
M67 observations. For each realization, we generate a 
population of synthetic single and binary stars (given an 
input binary frequency), and produce synthetic RVs for 
these stars on the true observing dates from our M67 
survey, which we analyze in the same manner as our real 
observations. We run this analysis for a range of true 
RV dispersion values, and for each value, we derive the 
difference between the synthetic combined RV dispersion 
and the input true RV dispersion. This difference is the 
co ntribution from un detected binaries (/? in Equation 5 
of iGeller et al.l[20^ . 

For this analysi s, we improve upon the technique of 
iGeller et al.l (j20ir)t) in our treatment of the binaries, and 
therefore, for clarity, we explain the method in some 
further detail here. In order to recover the true ve¬ 
locity dispersion, we must also estimate the binary fre¬ 
quency of the cluster. We use a s imilar Monte Carlo 
method to IGeller fc MathierJ (|2012ll to account for our 
incompleteness in bin ary detections ou t to th e hard- 
soft boundary. As in IGeller fc MathieiJ (j2012l) , we as¬ 
sume that the orbital parameters of the M67 binaries 
follow the same distributions as observed for the solar- 


type b inaries in the Galactic field from iRaghavan et al.l 
(|20ir)l) . ex cept here we use a circul arization period of 
12.1 days (iMeibom fc Mathiei^ l2005ll . We also i nake a 
few updates to the method of Geller fc MathieuI (|2012f) . 
First, rather than choosing one specific primary mass, 
we draw from an inferred primary-mass distribution 
within ou r M67 primary sa mple (derived by compar¬ 
isons to a lMarigo et al.l [200^ isochronel. Second, we at¬ 
tempt to recreate the correlation between binary mass 
ratio and orbital peri od observed for field binaries by 
IRaghavan et al.l (l2010l see Figure 17 and related discus¬ 
sion). Specifically, for binaries with periods between 100 
days and 1000 years, we choose mass ratios from a uni¬ 
form distribution with 10% of systems having a mass 
ratio of unity (e.g., twins). For binaries with periods 
greater than 1000 years we draw mass ratios from a uni¬ 
form distribution limited to be less than 0.95. For bina¬ 
ries with periods less than 100 days, we enforce ~19% 
to have mass ratios between 0.2 and 0.45, ^44% to have 
mass ratios between 0.45 and 0.9, and the remainder to 
have mass ratios between 0.9 and 1, all drawn from uni¬ 
form distributions between the respective mass ratio lim¬ 
its. Finally, we limit the field log-normal orbital period 
distribution so th at the binaries are d etached (using radii 
estimates from a lMarigo et al.l[20081 isochrone), and the 
periods are less than the hard-soft boundary. 

The assumed location of the hard-soft boundary, which 
is at roughly 10^-10® days in M67, is especially important 
for the undetected binary correction. We do not expect 
to detect binaries at these long periods. Furthermore, 
the hard-soft boundary is near the peak of the log-normal 
distribution, and therefore the assumed binary frequency 
is particularly sensitive to this cutoff. We estimate the 
hard-soft boundary as the location where a synthetic bi¬ 
nary’s binding energy is equal to the kinetic energy of a 
“typical” star moving at an assumed velocity dispersion. 
For the mass of this typical star, we take the mean mass 
of an object ( single or total mass of a binary) from the 
iHurlev et al.l (1200511 V-body simulation of M67, which 
at 4 Gyr is 0.95 Mq. We then assume a velocity dis¬ 
persion, run our analysis to derive the true velocity dis¬ 
persion in our M67 sample (given the resulting hard-soft 
boundary and total binary frequency), and iterate until 
this assumed velocity dispersion is < 0.05 km s“^ differ¬ 
ent from the derived true velocity dispersion in M67 (a 
somewhat arbitrary limit meant to be roughly equivalent 
to the precision with which we can measure the velocity 
dispersion). As a starting guess, we use the combined 
dispersion value of 0.80 km s“^, found above. 

Our analysis of the single main-sequence, giant and 
sub-subgiant stars in this sample requires only two itera¬ 
tions. The resulting orbital period distribution matches 
closely to those predict e d by the A^-body open cluste r 
models of iHurlev et al.l (j2(105ll and IGeller et al.l (1201,811 , 
and we estimate the total binary frequency in our sample 
to be 57% ± 4%. We derive a true velocity dispersion, 
after correcting for measurement error and undetected 
binaries, for the M67 single main-sequence stars, sub¬ 
giants and giants with V < 15.5 of 0.59jlg;Qg km s“^. 

Using a s ubsample of 20 of the brightest stars in 
our sample, iMathied (jl98,8ll measured a cluster RV 
dispersion, after correcting for measurement errors, of 
0.48 ± 0.09 km s“^. After also correcting for u ndetected 
binaries (assuming a 50% binary frequency), IMathied 
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(fl^ find an RV dispersion of 0.25 ± 0.18 km s 
I Girard et al.l (jl989ll show that this RV dispersion in¬ 
creases to 0.48±0.15 km s“^ (corrected for both measure¬ 
ment errors and un detected binare i s) wh en including the 
larger sample from iMathieu et id] (|1986f l. in good agree¬ 
ment with the r esult we find he re. _ 

More recently, iPasauini et al.l (|2012ll derive an RV dis¬ 
persion for main-sequence stars in M67 with V < 15 of 
0.680 ± 0.063 km s“^, and for M67 giants of 0.540 ± 
0.090 km s“^, using HARPS spectra with a typical pre¬ 
cision of '^lO m s“^ (not corrected for undetected bina¬ 
ries). If we divide our sample into the SM main-sequence 
stars with V < 15 and SM giant stars, after correcting 
for measurement errors, we find a combined velocity dis¬ 
persion for the main-sequence stars of 0.83 ±0.04 km s“^ 
and for giants of 0.60 ±0.08 km s“^. However, the main- 
sequence stars have a much higher binary frequency than 
the giants, and therefore a much larger correction for 
undetected binaries. After correcting for both measure¬ 
ment errors and undetected binaries, we find a true ve¬ 
locity dispersion for main-sequence stars (V < 15) of 
O.OOIq q® km s“^ and for giants of 0.60 ± 0.10 km s“^. 
Thus we find essentially identical velocity dispersions for 
the main-sequence and giant stars. We note that Padova 
isochrones indicate a mean mass for the main-sequence 
stars in this sample of 1.11 Mq (with a standard devia¬ 
tion of 0.10 M0), and for the giants of 1.33 Mq (with a 
standard deviation of 0.01 M0); our sample covers only 
a narrow mass range. 

One-dimensional proper-motion dispersion measure¬ 
ments from the literature are somewhat higher than 
these R V dispersions (as was als o note d bv iGirard et al.l 
Il989ll . iMcNamara fc Sanders! (|1978ll find a velocity 
dispersion of 0. 95 km s~^ with a Icr upper limit of 
<1.48 km s“^. IZhao et ^ (Il996ll find a d ispersion of 
0.96 ± 0.09 km s~^. and iGirard et al.l (|1989ll find a dis- 
persion of 0.81 ± 0 .10 km s“^. As discussed above, the 
I Girard et al.l (jl989l) proper motions have the highest pre¬ 
cision, and therefore here we will compare directly to 
their result. After correcting for undetected binaries, 
our true RV dispe rsion measurement differs from the 
I Girard et al.l ()1989[ 1 proper-motion dispersion measure¬ 
ment at about the 2a level (without accounting for an 
additional uncertainty on the proper-motion dispersion 
from the range in cluster distances quoted in the liter¬ 
ature). Given the uncertainties on these measurements, 
the different distances assumed in converting from angu¬ 
lar proper motions to km s“^, and the different methods 
used to account for measurement error in the proper- 
motion dispersions, we conclude that the RV and proper- 
motion dispersions are in agreement at our level of pre¬ 
cision. 

We have also examined the true RV dispersion of M67 
as a function of radius from the cluster center, shown 
in Figure [TOl We divide this sample of main-sequence, 
subgiant and giant single members into five equal bins 
in radius. Crosses, open circles and filled circles show 
the observed, combined and true RV dispersion values, 
respectively, at each bin in radius. The correction for 
undetected binaries is highest in the inner-most bin, as 
this bin has the highest binary frequency, consistent with 
our finding that the binaries are mass segregated with 
respect to the single stars (see Figure [9]). 



R (arcmin) 

Fig. 10. — Radial-velocity dispersion as a function of ra¬ 
dius from the cluster center. Only single main-sequence stars, 
subgiants and giants with V < 15.5 are used for this analy¬ 
sis. The observed radial-velocity dispersions are shown with 
crosses, combined radial-velocity dispersions (after correcting 
for the contribution from measurement error) are shown in 
open circles, and the true radial-velocity dispersions (after 
also correcting for undetected binaries) are shown in filled 
circles. The horizontal bars show the range in radius of each 
bin, and the vertical bars show the uncertainties on the ve- 
locity dispersion v alues. All uncertainties are derived as in 
iGeller et al.l (120101 1. 

The radial distribution of the true velocity dispersion 
within the parameter space covered by ou r survey is 
consis tent with an isothermal distribution. IZhao et al.l 
(|1996ll find an increase in the proper-motion dispersion 
as a function of radius, although they note that this effect 
may be in part due to increased field-star contamination 
with radius in their sample. We do not see a similar trend 
in our data. We have also examined different binnings 
with no detectable effect on the results presented here. 

Given the velocity dispersion, we can estimate the 
mass of the cluster through the virial theorem. Under 
the assumption that the cluster is in dynamical equilib¬ 
rium, the kinetic and potential ene rgies are related by 
1/2 _ 'qGM/R. From (jSpitzed[19871) . if this equation is 
rewritten in terms of the half-mass radius (ch), rj ^ 0.4 
for most systems. Also, for an isotropic velocity disper¬ 
sion, the observed half-mass radius in projection is given 
by rhp ~ and the one-dimensional RV dispersion 
is related to the rms velocity V by cr^ = U^/3. Thus 
the virial mass of the cluster can be calculated by the 
equation: 

Mv = 10 ^^. ( 2 ) 

Cr 

iFan et ahl (|1996D find a half-mass radius of M67 for 
stars with 13.8 < U < 14.5 of 2.54 ± 0.41 pc and for 
stars with 15.6 < U < 14.5 of 2.74 ±0.29 pc, assuming a 
distance of 850 pc. The M67 sample used here extends 
from 1/ 8 to U = 15.5, though the vast majority of 

the stars have V > 12.5. We will assume a projected 
half-mass radius for the stars in our sample of r^p = 
2.58 ± 0.45 pc. This value lies in the mi ddle and extends 
to the limits of the range in rhp found bv iFan et al.l (I1996D 
for these two magnitude regimes. Given our true RV 
dispersion above of ar = 0.59jlQ;gg km s“^, we find a 
virial mass for M67 of 2100l55g ^0- 

This value is in good agreement with previous dynam- 
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ical mass estimates. iMcNamara fc Sanders! (119781) de¬ 
rive a virial mass of 1600 Mq (with a Icr upper limit of 
4000 Mq), and a direct clu s ter m ass of 1100 Mq based 
on star counts. IZhao et (|1996H find a virial mass of 
1500 ± 250 Mq using their proper motions. From ex¬ 
amina tion of the cluster luminosity function, iFan et ;d] 
(|1996H calculate a mass of 1270 Mq within a radius of 
2000 arcsec and including stars with masses >0.5 Mq. 

These mass estimates (including o urs) are somewhat 
higher than the mass esti mates of Montgomery et al.l 
(1 19931) . iFrancid (|1989D and iMathieu ( 19831) . who find 
cluster masses of 724 Mq, 553 Mq and 903 Mq, re¬ 
spectively, all through analyses of the cluster luminosity 
function. We interpret these measurements as lower lim¬ 
its on the true cluster mass, as each of these surveys were 
limited either in magnitude or radial extent. 


7. SUMMARY 

This is the first in a series of papers studying the dy¬ 
namical state of the old open cluster M67 through precise 
RVs. Here we focus on determining cluster membership 
and identifying binaries (Sections [5]) within a complete 
sample of cluster stars. In total, we present results from 
13776 RV measurements of 1278 stars in the direction 
of M67 spanning from the brightest stars in the cluster 
down to about 4 magnitudes below the main-sequence 
turnoff (V = 16.5), covering a mass range of about 
1.34 Mq to 0.76 Mq, and extending spatially to 30 ar- 
cmin in radius from the cluster center (about 6 or 7 core 
radii). We combine RV measurements from multiple tele¬ 
scopes using different instruments, carefully accounting 
for their different single-measurement RV precisions (see 
Section IQ) . The vast majority of the stars in our sam¬ 
ple have multiple epochs of RV measurements, allowing 
for the detection of binary companions. For some stars, 
the time span of our observations exceeds 40 years, and 
about one third of these stars have observations spanning 
at least 10 years. 

In our analysis, we require at least 3 RV measure¬ 
ments before attempting to determine RV variability 
or RV membership probabilities. We also uti li ze th e 
proper-motion membership surveys o f ISandersI (Il977l). 
iGirard et al.l (|1989D . IZhao et all (| 19931) and lYadav et al.l 
( 20081 ) to remove additional non-members from our clus¬ 
ter sample. Using both proper-motion and RV member¬ 
ship information, we identify 562 cluster members within 
our sample, 142 of which show signihcant RV variability, 
indicative of a binary companion (or multiple compan¬ 
ions). 

We define a primary sample of stars with V < 15.5 
and within 30 arcmin from the cluster center, where we 
have >3 RV measurements for all proper-motion mem¬ 
bers (and, considering the entire primary sample, at least 
one RV measurement for all but 1 star and >3 RV mea¬ 
surements for all but 4 stars; see Section |4T]) . We then 
use this primary sample to construct a CMD cleaned 
from field-star contamination (Section |6T|), identify and 
discuss a few notable stars ('Section I6.2L compare the 
projected radial distributions of different cluster popula¬ 
tions (Section lOl) . and determine the RV dispersion and 
virial mass of the cluster (Section 16.411 . 

Within our cluster member sample, we identify 14 BSS, 
4 bright “yellow giants” residing on the GMD between 
the typical BSS region and the giant sequence (though 


one is outside of a 30 arcminute radius from the cluster 
center), and 2 “sub-subgiants” located to the red of the 
main-sequence but fainter than the subgiant and giant 
branches. These “exotic” stars have a remarkably high 
binary frequency of (at least) 80%, as compared to 22.7% 
detected binaries amongst the rest of the sample. 

Both the binaries and BSS are significantly more cen¬ 
trally concentrated than the single stars (i.e., non-RV 
variables) in our sample. Within our sample the bina¬ 
ries are, on average, more massive than the single stars, 
and therefore this result confirms that M67 is mass segre¬ 
gated (as also found by other authors), which is expected 
for a clust er that has lived th r ough tens of relaxation 
times. As IMathieu fc: Lathan] (|1986D also discuss, the 
central concentration of the BSS suggests that they too 
are more massive than the single stars in our sample, 
which is consistent with the predictions of theoretical 
formation channels. We do not observe a biinodal radial 
distribution for the BSS, as is observed for BSS in many 
globular clusters and those in the old open cluster NGG 
188. However, our sample only extends to roughly 6 or 
7 core radii, and a larger radial extent would likely be 
required to detect a possible halo BSS population and 
bimodal radial distribution. 

Finally, we determine the RV dispersion of the single 
main-sequence, subgiant and giant members in our pri¬ 
mary sample. Accounting for measurement errors, we 
find a combined RV dispersion of 0.80 ± 0.04 km s“^. 
When also corrected for the contribution from unde¬ 
tected binaries, we find a true RV dispersion of the cluster 
of 0.59jlo;o6 km s“^. The radial distribution of the true 
RV dispersion within our sample is consistent with an 
isothermal distribution. Using this true RV dispersion 
and a projected half-mass radius of Chp = 2.58 ± 0.45 pc, 
we calculate a virial mass for M67 of 2100jl55o ^0- 

Our long-term RV survey of M67 enables a detailed 
study of the cluster’s binary population, from the main- 
sequence through the giant branch and including the 
rich population of M67 BSS. Future papers in this se¬ 
ries will focus on the binary properties of the cluster 
(e.g., binary frequency and distributions of orbital ele¬ 
ments). In particular, we have determined orbital so¬ 
lutions for 108 binary members of the cluster, which 
we will present and analyze in subsequent papers. Ad¬ 
ditionally, we will study the binary properties of the 
BSS in detail, which are critical for our understanding 
of their formation mechanism (s) (e.g. iGeller fc MathieuI 
l20111 : IMathieu fc GelleHl20'T^ . Indeed, binary stars play 
a primary role in the dynamical evolution of star clus¬ 
ters and the formation of exotic stars like BSS. With the 
addition of the binary properties to the large body of 
existing observational work on the cluster, M67 will be 
invaluable to our theoretical understanding of stellar dy¬ 
namics, stellar evolution, the formation of BSS, and the 
long-term evolution of star clusters. 
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TABLE 4 

Radial-Velocity Data Table 
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RA 
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eji 
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Comment 
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8:51;22.91 

11:48:49.4 

12.720 
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0 

34 

45784.84 

47519.82 

33.32 

0.17 


98 

99 


99 

92 



BM 

SB2,CX111,PV 

2001 

S1027 

8:51:24.95 

11:49:00.8 

13.240 

0.600 

14 

13 

46808.93 

56707.92 
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0.15 

0.67 

93 

100 

93 

99 

95 

1.18 

0.000 

SM 
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13.260 
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2 

32 
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98 
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94 

99 

91 
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BM 

SBl 
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11:49:02.5 

15.210 

0.790 

17 

7 
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97 
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72 
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BM 

SB1,PV 

5001 

M5754 

8:51:23.50 

11:49:05.8 

16.220 

1.000 

3 

0 

56054.68 

56669.01 

34.24 

0.27 

0.74 

98 

96 




0.63 

0.669 

SM 


1002 

S1023 

8:51:26.84 

11:48:40.5 

10.540 

0.570 

0 

7 

41072.69 

50822.92 

3.92 

0.38 

0.49 

0 

84 

0 

7 

0 

2.04 

0.000 

SN 


2002 

S1040 

8:51:23.77 

11:49:49.3 

11.520 

0.870 

1 

57 

41073.67 

54164.75 

33.01 

0.08 

0.56 

98 

97 


99 

95 

10.05 

0.000 

BM 

SB1,X10,CX6,PV,YG 

3002 

S1018 

8:51:24.09 

11:48:21.9 

12.830 

0.570 

0 

32 

46874.69 

48291.90 

33.38 

0.16 


98 

100 

96 

99 

89 



BM 

SB2 

4002 

S1030 

8:51:25.95 

11:49:08.9 

13.230 

0.570 

2 

5 

46808.94 

56051.71 

34.22 

0.22 

0.77 

98 

100 


98 

90 

0.76 

0.653 

SM 


5002 

S1032 

8:51:26.52 

11:49:20.3 

13.480 

0.570 

2 

4 

47489.97 

55934.05 

34.24 

0.35 

0.75 

98 

100 


98 

96 

1.15 

0.169 

SM 

PV 


The contents of each column are defined in Section 5. 
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